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Abstract 

Self-interacting dark matter (SIDM) was introduced by Spergel & Steinhardt to address possible 
discrepancies between collisionless dark matter simulations and observations on scales of less than 
1 Mpc. We examine the case in which dark matter particles not only have strong self-interactions 
but also have strong interactions with baryons. The presence of such interactions will have direct 
implications for nuclear and particle astrophysics. Among these are a change in the predicted 
abundances from big bang nucleosynthesis (BBN) and the flux of 7-rays produced by the decay 
of neutral pions which originate in collisions between dark matter and Galactic cosmic rays (CR). 
From these effects we constrain the strength of the baryon-dark matter interactions through the 
ratio of baryon - dark matter interaction cross section to dark matter mass, s. We find that BBN 
places a weak upper limit to this ratio <J 10 8 cm 2 g _1 . CR-SIDM interactions, however, limit the 
possible DM-baryon cross section to < 5 x 10~ 3 cm 2 g _1 ; this rules out an energy-independent 
interaction, but not one which falls with center-of-mass velocity as s oc \ jv or steeper. 
PACS Numbers: 26.35.+c; 95.30.Cq; 95.35.+d; 98.70.Sa; 13.85.Tp 



1 Introduction 



Compelling observational evidence suggest that matter in the universe is dominated by "cold dark 
matter" (CDM), the simplest model of which is one where dark matter particles have interaction 
strengths at or below the weak scale, and thus interact today only through gravity (collisionless 
dark matter). The collisionless cold dark matter model has been proven through analytic and 
numeric simulations to be very successful in explaining the large scale structure of the universe Q . 

However, a series of observations of dark matter structures on the order of < 1 Mpc (e.g., 
galaxy rotation curves & [|, ||], strong lensing from galaxy clusters ||, Tully-Fisher relation of 
spiral galaxies bar stability in spirals 0, |8| and deficit of low-mass subhalos in the Local Group 
|]]) appear to contradict the prediction of collisionless CDM simulations that dark matter halos 
form high-density cores, favoring shallower, lower-density dark matter halo cores. 

Spergel and Steinhardt |[(]] have proposed a variation of CDM in order to alleviate these dis- 
crepancies. In their proposed picture, dark matter particles interact with each other strongly and 
have a large scattering cross section, but negligible annihilation or dissipation.^] Subsequent numer- 
ical simulations |11, ^] have shown that this collisional cold dark matter, or SIDM, does indeed 



predict halo cores in better agreement with the observations mentioned above, provided 

0.5 < s < 6cm 2 g -1 . (1) 

where the the parameter s = 0"sidm/-^d is the dark matter-dark matter elastic scattering cross 
section over the dark matter particle mass. Potential difficulties for this model may arise at cluster 
scales Jl3|, [|, IHf ; one way to address these would be a velocity-dependent cross section, the 
effect of which we will consider below. 

If dark matter particles interact with each other through the strong force, then similar inter- 
actions would be expected between them and ordinary baryons as well, with cross sections of the 
same order ]12|] . Wandelt et al. considered this case and found, perhaps surprisingly, that such 
interactions cannot be excluded by galaxy halo data nor from direct observation by space-borne 
cosmic-ray detectors, if dark matter particles have a mass larger than 10 5 GeV. Earlier papers ]i~6[ 
have also considered astrophysical constraints on SIDM both with and without interactions with 
baryons. 

From the particle physics point of view, we choose a model-independent approach. The ques- 
tion we address is "If a massive neutral particle exists, how can we constrain its properties from 

1 If the dark matter particle and antiparticle states were equally populated and annihilations were allowed, these 
would proceed vigorously, leaving a relic abundance today Q,n ~ 10 -37 cm 2 {cr ann u/c) _1 which would be negligible 
for the case of strong interactions. This implies that the SIDM particles either have strongly suppressed annihilation 
cross sections, or are asymmetrically populated (as is the case for, e.g., baryons). 
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astrophysical observations?". We refrain from a detailed discussion of the naturalness of such a 
particle except to note how quickly model builders proposed several particle physics implementa- 
tions (e.g. ||17|| ) of self- interacting dark matter. Some promising candidates, the 50 and low-mass 
strangelets, were then ruled out j^]. Nevertheless, there remain intriguing candidates for dark 
matter which lie in the regions of masses and baryon/self- interaction cross sections we consider, 



both within super symmetric extensions of the standard model [18] or as possible relics from the 



QCD phase transition within the standard model [19]. Our results place constraints on the energy 
dependences of the cross sections in these models. 

In this paper, we examine the effect of inelastic SIDM-baryon interactions. Specifically, we 
consider the impact of such interactions on big bang nucleosynthesis (BBN), and on the 7— rays 
produced by interactions of cosmic ray nuclei with dark matter particles. In BBN, the formation of 
the light elements deuterium, 3 He, 4 He, and 7 Li, depends on the time at which the weak interactions 
"freeze out," the free neutron decay rate, and the temperature at which nucleosynthesis occurs, 
the latter being determined by the photo-dissociation of deuterium. If dark matter can interact 
strongly with baryons, it is possible that it can destroy light nuclei. Since deuterium has the 
smallest binding energy and determines the onset of nucleosynthesis, any dark matter deuterium 
dissociation will change where the onset of nucleosynthesis happens and thus provide a sensitive 
probe to baryon-dark matter interactions during BBN. We will see that the constraints provided 
by BBN are quite weak-in other words, for the parameter ranges of interest, SIDM has negligible 
effect on BBN despite the strength of the SIDM-baryon interaction. 

The interaction of cosmic rays with interstellar matter leads to observable signatures in the 
form of 7— rays. Collisions between energetic cosmic ray protons and interstellar hydrogen can 
produce neutral pions, which in turn decay into two high-energy 7— rays; this emission is believed 
to dominate the observed flux of high-energy photons from the Galactic disk 2C, 22 1. Including 
dark matter-baryon interactions allows for similar processes to occur between cosmic ray nuclei 
and ambient SIDM, thus possibly contributing to the 7— ray diffuse background. The 7— ray flux 
from CR-SIDM interactions will be a function of the assumed cross section for the interaction 
between dark and baryonic matter and can be calculated given the functional form of the Galactic 
cosmic ray flux and dark matter distribution. The result can then be compared with diffuse 7-ray 
observations from 7— ray telescopes (most recent of which are the EGRET observations, [21, p0| ) 
and strong upper limits can be placed on the cross sections. 

This paper is organized as follows. We will discuss the strong interacting dark matter cross 
section in §|2[ In §^ we determine the SIDM effect on BBN and its constraints. In §|] we compute 
the expected 7-ray flux from cosmic ray - SIDM interaction in the Milky Way and the resulting 
constraints on the cross section. Other possible constraints are noted in §[|. Finally, our findings 
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are summarized and discussed in §pl 

2 SIDM-Nucleon Interactions 

We wish to characterize the inelastic interactions between baryonic nuclei and SIDM particles D. 
We will express the cross sections for inelastic interactions in terms of the elastic cross sections. 
The inelastic behavior is assumed to obey simple scalings with the masses of the reactants flO| 
and be energy-independent (but we will use the different constraints below to infer constraints on 
possible energy dependence of SIDM-baryon interactions). 

We begin by considering the SIDM-nucleon interaction, with an energy-independent elastic 
scattering cross section a^ff 10 . This need not be the same as the SIDM self-interaction cross 
section, but was taken to be so in the work of |l~2|| . One can then generalize this to determine the 
cross section for elastic SIDM- nucleus scattering. Following 16], we adopt the form 



.elastic A 2 { K A )\ .elastic /r>\ 

°DA =A (^yj a DN (2) 

where A is the atomic mass number of the interacting nuclide, and fJ,(A) is the reduced nucleus-dark 
matter mass. Finally, we will be interested in inelastic scattering, i.e., when internal degrees of 
freedom are excited in the baryonic component. We will write the (energy-independent) inelastic 
scattering cross section for channel % as 

i „ .elastic /n\ 

where Oj encodes the relative strength of the inelastic interaction. 

A typical strong interaction has a cross section, £j|^y tlc ~ 1 barn and it is assumed that a dark 
matter particle has a mass, Mb <^ 1 GeV. A useful parameter given these typical values for the 
cross section and mass is defined as 

s = atf ic /M D (4) 

usually quoted in cgs units as cm 2 g _1 . A useful conversion factor for particle physics units is 
1 cm 2 g- 1 = 1.78 barn/(GeV/c 2 ). 

3 Limiting SIDM with BBN 

Primordial nucleosynthesis is the process by which the light elements, consisting mainly of deu- 
terium, 3 He, 4 He, and 7 Li, were produced. Standard BBN is a one parameter theory, depending 
only on the baryon-to-photon ratio, rj = ii-q/ii^ = 2.74x 10 _8 J1b^ 2 > where ^b is the current baryon 
density relative to the critical density, p crit = 3H$/8ttG, and H = lOOh km s _1 Mpc -1 is the 
present value of the Hubble parameter. 
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Before considering the effect of SIDM, it is useful to recall the basic physics of BBN. At the 
onset (T > 1 MeV; t < 1 s) of the BBN epoch, the universe was dominated by relativistic particle 
species, initially consisting of photons, 3 species of light neutrinos, and pairs; baryons were 
non-relativistic and essentially only n and p with no complex nuclei. At this stage, all of these 
particle species existed in thermal and chemical equilibrium, due to the interaction rates between 
particles being much faster than the expansion rate of the universe. 

At a temperature around an MeV, and time about 1 s, the weak interactions, coupling the 
neutrinos to the other particle species and keeping the n and p in equilibrium, become ineffective 
(interaction rates less than H) . At this time the neutrinos fully decouple from the primordial plasma 
and the neutron abundance "freezes out", relative to the proton abundance, with the exception 
of the occasional free neutron decay. Shortly after, the electrons and positrons start to become 
non-relativistic, and thus annihilate, producing photons, which in turn heats the plasma relative 
to the neutrinos. Trace amounts of the light nuclei, deuterium, 3 He, 4 He, and Li exist, with 
abundances held below nuclear statistical equilibrium levels due to deuterium photo-dissociation 
by photons in the high-energy tail of the thermal distribution. This "deuterium bottleneck" persists 
until the universe cools to T ~ 0.07 MeV, when energetic photons are sufficiently rare, at which 
point deuterium rapidly grows and then is burnt into the other light elements. 

We now consider what role CDM plays in these series of events. Depending on the origin of the 
CDM and its coupling to Standard Model particles, CDM may or may not be in thermal equilibrium 
with the rest of the cosmic plasma. In the standard CDM scenario, the CDM only interacts weakly 
with itself and with Standard Model particles, and thus these interactions will freeze out well before 
the epoch of BBN. Then the dark matter cools rapidly due to Hubble expansion. Throughout BBN, 
the dark matter is non-relativistic and the Universe is radiation dominated. So in the standard 
case, dark matter plays no role at all. 

What happens in the SIDM case? The dark matter is still non-relativistic, but now elastic 

collisions between baryons and dark matter will thermalize the dark matter at early times and 

bring Tb = T 7 . Dark matter-baryon decoupling will not occur until much later, and we note 

that this will potentially have interesting effects on structure formation^. However, for our study 

2 The tight coupling of the baryon and photon fluids and the coupling between the dark matter and the baryons 
result in coupling the DM to the photons. This gives pressure to the dark matter fluid and changes the evolution 
of perturbations until the dark matter decouples from the baryons and then cools adiabatically. We conjecture 
that this effect leads to a suppression of fluctuation power on perturbation modes which entered the horizon before 
DM-baryon decoupling. Standard freeze-out arguments show that DM-baryon decoupling occurs at a redshift of 
z~2x 10 4 <r^f ic _(3/2) , where 

^.clastic j g measurec j m barn. The characteristic mass scale of perturbations entering 
the horizon at this time is of order 10 9 Mq(u% !1 1 !} Uc ) 2 . Note that this scaling now depends just on the cross section, not 
on the mass of the particle. This suggests that DM-baryon interactions reduce the dark matter perturbation power 
on galaxy scales — just those which motivated the study of interacting dark matter in the first place. We plan to 
study these effects in more detail in a future publication. 
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of BBN elastic scattering is not important except for thermalizing the dark matter. What does 
affect BBN are inelastic interactions with SIDM particles. Specifically, such collisions can impart 
internal energy to the nucleus and lead to its dissociation. Given that deuterium has the smallest 
binding energy and its photo-dissociation determines the onset of nucleosynthesis, we will focus on 
deuterium breakup; similar processes should occur for other nuclei but are not significant due to 
the higher binding energies. 

We consider the reaction D+d — > D+n+p and neglect the rare, three-body, inverse reaction. In 
order for this interaction to happen the center-of-mass kinetic energy of the dark matter-deuteron 
system must be equivalent to at least the binding energy of the deuteron, B = 2.224 MeV. We 
also assume that the break-up cross section is comparable to the scattering cross section, a D k f = 
a i a D<f tlC ' where ot\ represents the efficiency of the break-up reaction and is of order unity. We 
relate the deuterium-dark matter scattering cross section to the proton-dark matter cross section 
as defined by eq. ). 

In making this interaction applicable for BBN, we need to compute the thermally averaged 
reaction rate, (crv), where v is the relative velocity of the interacting particles. Since for our 
discussion, we consider a constant cross section with the appropriate energy threshold, this thermal 
average can be easily found. 

M = 5.221 x lOVaf^ { ^ 3/2 (l + f ) Tl' 2 exp (-^) cmV 1 (5) 

V 1 + tie) 

where m p and rrid are the masses of the proton and deuteron respectively, and Bg and Tg are the 
binding energy of deuterium and the temperature expressed in units of 10 9 K.0 

The addition of this dark matter reaction adds another channel for deuterium destruction in 
addition to photo-disintegration. Thus, it will take longer to build up enough deuterium so that 
nucleosynthesis can occur. To study the effects of the new interaction analytically, we will look at 
the quasi-static equilibrium (QSE) abundance of deuterium, as laid out by, e.g., |23|, 24], and see 
how this changes when nucleosynthesis occurs. The abundance of deuterium will be most sensitive 
to the np —¥ dj, dj — > np and dD — > Dnp reactions. The QSE abundance is then given by the 
ratio of sources to sinks. 

v.QSE = „ /mj = ^jkiYjYkjjk ^ il] Y n Y p [np -> dj\ 

1 ~ l/ B E jk i Yi [U - 3 k] d Y^ [d 7 -> np] + Y D [dD - Dnp] { > 

Here rtj is the i th particle species number density, ne is the baryon number density, Yy = I/77, and 

3 Throughout this analysis we assume that the dark matter remains non-relativistic during the epoch of BBN, 
Md 10 MeV. If dark matter is relativistic during this time, not only will the reaction rate take on a new form, but 
now the dark matter will contribute to the expansion of the universe. 
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[jk — ► il] are the reaction rates 

[jk -> il] = N A pB{<Tjk->u v) (7) 

given by several compilations, where Na. = fn~ l is Avogadro's constant, and ps is the baryon mass 
density. 

The forward and reverse reaction rates are related to each other by detailed balance. This tells 
us that the ratio of reverse to forward rates is given by the thermal equilibrium distributions of the 
interacting species. For np •<=>■ 7c? one finds the Saha expression 

1^^^ y^EQ^EQ 5 -3/2 , r> irp \ / R x 

[ np _> dl ] ~ yJQyEQ " L4 ° X 10 T » exp ^~ B ^ T ^ ^ 
Assuming dark matter has a mass Mjj ^> m p , the rate is given by the following expression. 

v- D[rfD . D „^, 221xl „v(^)(^)(-^)( 1 + |)rr»p(-f ) (9) 

The ratio of the two deuterium destruction rates is then 



Y 7 [dj — > np] V cm 2 g 1 / \ . 

assuming Tg <C Bg and where r/10 = r//10 — 10 . For reasonable values of a\ and s, we see that this 
ratio is tiny for the temperatures relevant to BBN. This simply reflects the smallness of 77, i.e., the 
large photon-to-baryon (and photon-to-SIDM) ratio. The smallness of the SIDM contribution will 
lead to a very weak perturbation to the light elements. 

Combining the results of eqs. d), (|), and © and requiring = gives the QSE deuterium 
abundance. The deuterium bottleneck ends, and the light elements are formed, when this increases 
to order unity, at the temperature 

25.81 



(11) 

34.18 - In (7710) - § In (Tg) + 2.1 x 10- 9 r/ 10 ai (s^t) (g&) Tg 
with Tg <i 1. 

Thus one can see that with ~ 10, one needs an extremely large cross section, a\s > 10 
cm 2 g _1 , to cause a noticeable change in the light element predictions. 4 He is the most sensitive 
to these changes. To derive a limit on SIDM, we will adopt a very generous observational lower 
bound Y p > 0.230 [p5f|. Using this, and a full numerical implementation of eq. @ in the expanded 
BBN code discussed in [p6fl , yields a limit a^s < 1.263 xlO 8 cm 2 g _1 for the case Mq 3> m p . Given 
the rate dependence on dark matter mass and cosmological parameters we can derive an explicit 
form for the maximum cross section, 

a lS < 1.263 x 10 9 ^ ( cm ' g _1 (12) 
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for a baryon-to-photon ratio r/io = 5.8 and requiring Y p > 0.230. The constraints are only slightly 
different for different assumed rj and helium abundances. 

We see, therefore, that BBN cannot place strong constraints on dark matter mass or strong 
baryon cross section. In other words, energy-independent SIDM is completely compatible with light 
element constraints and BBN. The weakness of the constraint follows from the smallness of r\ and 
its imposition of the deuterium bottleneck which delays nuclear reactions until temperatures far 
below the binding energies of the light nuclei. This ensures that SIDM-deuteron interactions will 
be too weak to allow dissociation, and the lack of SIDM-nucleon bound states ensures that SIDM 
particles will not otherwise compete with nucleons as the light elements are built up. 

Having seen that SIDM-nucleon interactions are compatible with nucleosynthesis in the early 
universe, we now turn to a present-day consequence of these interactions, namely the production 
of 7-rays from SIDM collisions with cosmic rays. 

4 Limiting SIDM with Cosmic Rays 

The Galaxy is penetrated by a flux of nonthermal, relativistic nuclei, the cosmic rays. The bulk 
of these particles (i.e., those with E < 100 TeV) are of Galactic origin, and have kinetic energies 
typically ~ 1 GeV with a power law energy spectrum 

where 7 ~ 1.7. Inelastic collisions of cosmic ray nuclei (primarily protons) with SIDM particles are 
detectable in principle by the gamma-ray signature of the decay of the neutral pions produced in 
such interactions: 



D + p-> D + A + 



P+ 7T° 



^ 7 + 7 

where D is the dark matter particle. This reaction is simply the SIDM analog of the usual pion 
production pp — > ppir — > 77 reaction, responsible for the bulk of the Galactic contribution to the 



7-ray sky 127L B8L 21 



The Galactic gamma-ray signature for this process is readily calculated. We assume axial 
symmetry for the cosmic-ray and SIDM spatial distributions, and a relevant dark matter-nucleon 
cross section <7^J- astlc independent of the incoming proton energy. The emissivity g 7 (7— ray photons 
per volume per time) is then 

qi (r,z) = 2n D (r,z)a i Bt stic ^ P (r,z) . (14) 

Here, the factor of 2 is the number of photons per 7r° decay, n D (r, z) is the number density distri- 
bution of the dark matter in the Galaxy, and <3? p is the angle-integrated cosmic ray proton intensity 
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for all energies above the pion production threshold, 

POO 

$ p = 47t/^ <f> E ,pdE (15) 

(16) 



with 

m 7T o(2m p + m n o) 



Emm = (m p + m^o) 



1 + 



2M D (m p + m n o) 

For Mb larger than ~ 1 GeV, the second term in the brackets of eq. (|I6|) is 1 and E m m does not 
change appreciably with Md, and we have E m in ~ m p + m n o. Using the parameter s = afffi^ /Mp, 
and the fact that with typical dark matter halo speeds v/c ~ 10~ 3 <C 1, then po = M^ud and we 



can re- write eq. (14) as 

q 7 (r, z) = 2 p D (r, z) a 2 s $ p (r, z) . (17) 

The parameter a 2 as used in eq. (|T7|), is the product of the ratio a Dp^DNw / 'cr^v* 10 and the branching 
fraction /o for the decay of the A + to p + tt° rather than n + 7r + . From isospin considerations it 
follows that /o = 2/3. In the absence of a detailed particle physics model for the interacting dark 
matter, it is natural to expect that the CR inelastic collision strength is related simply to the 
breakup collision strength by just this isospin factor, a 2 = 2/3 a±, as both reactions are simply 
exciting internal degrees of freedom in the proton. 

The resulting 7— ray intensity 4>-y, DN in any given direction will then be the line integral of g 7 
along the line of sight: 



1 

4-7T Ao.s. 



b ~/,DN = T~ / 1l( r ' Z ) M ( 18 ) 



which, along a line of sight defined by a set of Galactic coordinates (£, b), takes the form 

(j) J>DN = — J°° g 7 ( ^J(R cos b cos £ - L) 2 + (R cos b sin I) 2 + (R sin b) 2 , R sin b) dR (19) 

where L is the distance of the Sun from the Galactic center, and R is the heliocentric distance of 
a point of Galactocentric distance r. 

A first estimate of the expected 7-ray intensity from cosmic ray - dark matter interactions can 
be made if we approximate the spatial distribution of the CR flux by a step-function in space: 
Uniform and equal to the demodulated solar value inside an oblate ellipsoid with a semi-major axis 
of 10 kpc in the Galactic plane and 1 (5) kpc towards the Galactic poles, and zero outside this 
ellipsoid. In the limit Md — > 00, eqs. (|l5|) and (16) give $ Pi0 o = H-8 cm" 2 s" 1 , while if we relax 



1 + 



m i o(2m p +m ir o) 



2M D (m p +m 7T a) 



our assumption of large Mjj, $ p will scale as <^ p = $ Pi0 o 

In addition, we assume that the mass distribution of the dark matter halo can be described as 
a flattened, non-singular isothermal sphere j29| ; p(r, z) = Pc^ 2 ./ (r^ + r 2 + (z/q) 2 ) , where p c is the 
central density, the core radius and q the flattening of the halo. In the following calculations we 
have used q = 0.8, p c = 0.05 M pc -3 and a value of = 5 kpc, such that the total dark matter 
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Direction of 
Observation 


SIDM expectation 

( ol^s \ photons 


EGRET observation 

photons 


implied «2 S 
upper limit 


\ 1cm 2 /g J cm 2 sr s 


cm 2 sr s 


Galactic center 


2.1 x ict 1 


5.5 x 10" 4 


2.6 x 1(T S 


Galactic anti-center 


6.7 x 10~ 3 


1.5 x 10" 4 


2.2 x 1(T 2 


Galactic poles 


2.5(12.) x 1(T 3 


1.3 x 10" 5 


5.2(1.1) x 10" 3 


M31 


9 x 10" 2 


< 8 x 10~ 5 


9 x 10" 4 



Table 1: Expected 7-ray intensities from SIDM-CR interactions and implied cross-section upper 
limits. The numbers in the parentheses indicate results for a semi-major axis of the cosmic ray halo 
ellipsoid in the direction perpendicular to the Galactic plane equal to 5 kpc, while other results 
refer to a semi-major axis of 1 kpc. 



mass included in a radius of 100 kpc is ~ 10 12 Mq. The distance of the Sun from the Galactic 
center was taken to be L ~ 8.5 kpc. 

We have calculated </> 7 in this simple approximation for directions of observation towards the 
Galactic center (£ = 0°, b = 0°), towards the anti-center (£ = 180°, 6 = 0°) and perpendicular to 
the Galactic plane (b = 0°) and we have compared the results with the corresponding EGRET 
observational results from Hunter et al and Sreekumar et al p^] to derive upper limits for 02 s. 
Our results for the limit of large Mjj are shown in Table 1. The general result for the cross section 
upper limit will be given by 



a2S < 02^00 



1 + 



m 7r o(2m p + m^o) 
2M D (m p + m^o) 



(20) 



where Sqo is the value as Md — * 00. 

Note that for a sufficiently flattened cosmic ray halo (1 kpc minor axis) the strongest constraint 
comes not from the direction where the 7— ray intensity is minimum (the Galactic poles), but from 
the observations towards the Galactic center, where the 7— ray intensity is maximum. 

For the derivation of the upper limits presented above, the entire 7— ray flux detected by EGRET 
was attributed to pion decay from cosmic ray-SIDM interactions. This is a good assumption for the 
direction towards the Galactic poles where the contribution of the 7— rays from collisions of cosmic 
rays and interstellar baryons is not dominant. However, in the Galactic plane most of the observed 
intensity can be attributed to the presence of interstellar baryonic gas. Thus, the constraints in 
the directions of the Galactic center and anti-center could be made even more stringent than these 
conservative limits if one introduces a model for the emission from the interstellar medium. 

We note that constraints of this type can be derived for external galaxies as well. In particular, 
we consider M31, which has the advantages of being nearby, and similar to the Milky Way. M31 
has been searched for in EGRET sky maps, but no positive detection has been possible and an 



upper limit of 1.6 x 10 cm s has been placed instead on its diffuse 7— ray flux 



Given 
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the measured dark mass Mp 31 = 6 x 10 10 M© in the inner 10 kpc of M31 pl| . we predict a SIDM 
7- ray flux of 

jMsl $M31 M M31 

C 1 = 2 " 4nR 2 ° (21) 

47m M31 



where i?M3i = 750 kpc is the distance to M31. Following [32|, we estimate the M31 cosmic ray flux 
<I>M3i by assuming supernovae to be the site of cosmic ray acceleration, and thus taking the cosmic 
ray flux to be proportional to the supernova rate. This gives <I>^ 131 = 0.453>^ w , which allows us to 
predict $^ 31 = 1.8 x 10 -5 a 2 s cm -2 s _1 (where s is given in cm 2 /g). This must be lower than the 
observed limit, which implies a limit 

a 2 s M31 < 9 x 10- 4 cm 2 g" 1 (22) 

Thus we see that the M31 7-ray limits are about an order of magnitude more constraining than 
those of the Milky Way; however, since they rely on more assumptions regarding the cosmic ray 
flux in M31, we regard the Milky Way limits as more secure. 

Finally, it is very intriguing to note that a recent reanalysis [^] of the EGRET data finds 
evidence for a diffuse Galactic 7-ray halo. Specifically, a wavelet analysis was used to identify a 
7-ray component at large angular scales, centered on the Galactic center, the intensity of which 
significantly exceeds the predictions of a model which includes known Galactic sources. One should 
bear in mind the error budget is large: the inferred halo intensity levels are uncertain to within 
a factor of ~ 2, and possible systematic effects remain. Also, it is entirely possible that the halo 
can be explained via inverse Compton scattering of background photons off of energetic cosmic ray 
electrons p2|. 

Nevertheless, it is tantalizing to interpret this diffuse halo in terms of SIDM-baryon inter- 
actions. The 100 MeV halo intensity levels toward the poles is quite uncertain but of order 
~ 10~ 6 cm~ 2 s _1 sr _1 , i.e., about 10% of the full EGRET result. In terms of our model, this 
would amount to a measurement of a 2 s at about 10% of our limit in Table 1, namely a 2 s ~ 
5(1) x 10 -4 cm 2 g . Future 7-ray observatories such as GLAST can verify the existence of this 
diffuse Galactic halo, and will be in a position to measure its spectrum. If indeed the flux arises 
from 7T° production in SIDM-baryon interactions, we would predict that the spectrum should reflect 
this origin, in particular, it should be a smooth continuum which shows the "pion bump" feature at 
E~ = m 7r o/2 = 67.5 MeV. In addition, similar halo emission should be visible around Local Group 
galaxies-not only M31, but also the Magellanic Clouds. 
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5 Other Possible Constraints 



SIDM-nucleon interactions can have other effects as well. The possibility of DM-baryon interactions 
contributing to the gamma ray background provokes the question whether some fraction of the ~ 20 
unidentified EGRET sources at high Galactic latitudes could be due to dark matter concentrations 
in the Milky Way halo. A rough estimate shows that typical gamma-ray fluxes of these sources 
would be compatible with DM clumps of order 10 3 Mq (R/l kpc) (1 cm 2 g _1 /a2s). While this is an 
intriguing speculation, the small number statistics of these sources makes distinguishing possible 
candidates from a background of approximately isotropically distributed AGN a difficult task. A 
more detailed study will therefore have to await the higher sensitivity and resolution data from 
future missions such as GLAST. 

Another potential (but model-dependent) constraint comes from diffuse radiation. Consider the 
effect of SIDM particle passage through the interstellar medium of our Galaxy, most of which is 
hydrogen. Elastic collisions between SIDM particles and the protons occur, with a relative velocity 
given by v ~ 200 km/s ~ 10 _3 c, and thus a center-of-mass energy Eqm — m p v 2 /2 = 0.2 keV. 
Given that we expect m p <C Md , it follows that a sizable fraction of this energy will be imparted to 
the recoiling proton. As a source of heating, this process does not violate observed constraints on 
the properties of molecular clouds [16|. However, if the scattering process leads to photon emission 
by the proton, this radiation might be observable. One can show that if each scattering event 
produces a photon with Eqm, then the resulting diffuse emission is at or above the level observed 
in soft X-rays. Of course, it need not be the case that this radiative scattering occurs, but given a 
particular SIDM model, this constraint should be investigated. 



6 Discussion and Conclusions 

We have considered constraints on strong interactions between dark matter particles and baryons. 
The big bang nucleosynthesis (BBN) constraints are weak, indicating that interacting dark matter, 
as it has been recently proposed, is fully consistent with BBN. However, very strong limits come from 
considering cosmic ray (CR) interactions with the dark matter. For dark matter masses Mp ^ 1 
GeV the ratio s = odnIMd is limited to be less than 3.9 x 10~ 3 cm 2 /g. Constraints from direct 
detection experiments such as the quantum calorimeter XQC |34|] require Md > 10 5 GeV(l2|. |35|| . 
Our limit therefore constrains dark-matter baryon interactions to be less than a hundredth of 
the self-interaction strength which, simulations suggest, is required to affect the structure of galaxy 
halos within the self-interacting dark matter (SIDM) scenario. Any particle physics model of SIDM 
must explain this order of magnitude suppression of the dark matter-baryon interaction. 

For the case in which the inelastic cross section is independent of energy, our results are sum- 
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Figure 1: Limits on the cross section-to-mass ratio s as a function of SIDM mass. The values 
needed for galaxy halos (eq. [l]) are shown in the shaded region. Excluded regions are hatched; 
single-hatched regions are excluded due to the inelastic interactions with Galactic cosmic rays; the 
cross-hatched regions are also excluded via BBN. The dashed curve shows the limit using M31 
7-ray observations. 
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Figure 2: Limits on the cross section-to-mass ratio s as a function of center of mass energy, for the 
case in which Mp ^> m p . The dashed curve shows the limit using M31 7-ray observations. 
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marized in Figure |]. It is of course possible, and indeed plausible, that SIDM-baryon interactions 
might be energy-dependent. Thus, we have displayed in Figure || the constraints on s as a function 
of the center-of-mass energy at which those constraints apply. Both figures show that energy- 
independent DM-baryon interactions are ruled out by the cosmic ray constraints, unless the ratio 
of DM-baryon scattering to DM self-scattering a < 10~ 3 . On the other hand, if the interaction 
is energy-dependent, Figure Q constrains the behavior. For example, if we assume that s oc Eq M , 
where Ecm is the center of mass kinetic energy, then we find using our limits from cosmic rays com- 
ing from the Galactic center that n < —0.39. Using the M31 constraint steepens this to n < —0.47. 
In both \jv scaling (or stronger dropoff) is allowed. Thus, our 7-ray constraints provide 

important information about the energy dependence of a putative SIDM-baryon interaction. If 
such an energy dependent cross-section were characteristic of dark matter self-scattering as well, 
SIDM would be less interactive in objects with high velocity dispersion such as clusters of galaxies, 
possibly softening constraints on SIDM from the statistics of strongly lensed arcs fl3|j . 

We stress that our results do not rely on any detailed assumptions about the underlying particle 
physics model governing the interactions. In the context of a more detailed particle physics model 
for SIDM, one could also examine signatures of less generic processes, for example the generation 
of diffuse emission of soft X-rays from radiative scattering between SIDM and interstellar medium 
particles. 

Finally, we note that elastic SIDM-baryon interactions thermalize the dark matter at early 
times. This may have potentially important consequences for structure formation. We plan to 
investigate this issue in more detail in a future publication. 
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